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1 Introduction 



One of the most important areas of high energy physics comprises tests of the theory of 
strong interactions, Quantum Chromodynamics (QCD) [1]. Since the theory contains 
in principle only one free parameter, the strong interaction scale A, tests of QCD can 
be quantified in terms of comparison of measurements of A in different processes 
and at different hard scales Q. In practice most QCD calculations of observables 
are performed using finite-order perturbation theory, and calculations beyond leading 
order depend on the renormalisation scheme employed, implying a scheme-dependent 
A. It is conventional to work in the modified minimal subtraction scheme [MS 
scheme) [2], and to use the strong interaction scale for five active quark fiavours. 
If one knows Aj^ one may calculate the strong coupling as{Q^) from the solution 
of the QCD renormalisation group equation [3]. Because of the large data samples 
taken in e"'"e~ annihilation at the Z° resonance, it has become conventional to 
use as a yardstick as(Af|-), where Mz is the mass of the Z° boson; Mz ~ 91.2 
GeV [4]. Tests of QCD can therefore be quantified in terms of the consistency of 
the values of as(Af|-) measured in different experiments; such measurements have 
been performed in e"'"e~ annihilation, hadron-hadron collisions, and deep-inelastic 
lepton-hadron scattering, covering a range of from roughly I to 10^ GeV^. In 
addition to testing QCD, the precise measurement of as{M'^) allows constraints on 
possible extensions to the Standard Model (SM) of elementary particles; see eg. [5]. 

Over the past decade many measurements of have been presented. These are 
reviewed here, with emphasis on the more recent results. The systematic uncertainties 
associated with these measurements, both experimental and theoretical in nature, are 
discussed, and the degree of consistency between as(Af|-) values is examined. The 
task of preparing this article has been aided by a number of previous reviews [3, 6]. 

2 Theoretical Considerations 

An inclusive observable X may be written schematically: 

X = (I + ,5^^^) (I) 

where X^^ represents the electroweak contribution. Since, with observables of this 
type, as enters via the small QCD radiative correction, S'^'^'^ , a precise measurement 
of as generally requires a large data sample. Observables can also be defined that are 



2 



directly proportional to S'^'^'^ and hence potentially more sensitive to a^. In either 
case can be separated into perturbative and non-perturbative contributions: 

^QCD ^pert _j_ ^non-pert 

The perturbative contribution can in principle be calculated as a power series in 
cts, though in practice the large number of Feynman diagrams involved renders 
calculations beyond the first few orders intractable. An observable must be calculated 
to at least next-to-leading order in the MS scheme so as to define Aj^; the solution 
of the renormalisation group equation [3] to the same order can then be used to 
translate consistently to the yardstick as{M^). The non-perturbative contribution, 
often called a 'hadronisation correction' in e"'"e~ annihilation or a 'higher twist effect' 
in lepton-hadron scattering, is expected [7] to have the form of a series of inverse 
powers of the physical scale. Hence 

= + SA. ,3) 

In general, due to the presence of 'renormalon ambiguities' in perturbation theory 
[7], the perturbative contribution cannot be calculated independently of consideration 
of the non-perturbative power-law contribution. For several inclusive observables 
the perturbative series has been calculated to next-to-next-to-leading order; for 
jet-like observables it has been calculated to next-to-leading order. The power-law 
corrections cannot in general be calculated, though there has been recent progress in 
this direction [7]. Therefore, in any comparison of a QCD prediction with data, the 
uncertainties relating to both the uncalculated higher-order perturbative, as well as 
non-perturbative, contributions should be estimated, and a theoretical uncertainty on 
the extracted value of as{M'^) assigned accordingly. 

From an operational point-of-view, truncation of the perturbative series at finite 
order causes a residual dependence on the scheme-dependent renormalisation scale 
fi. This parameter is formally unphysical and should not enter at all into an 
exact infinite-order calculation, and its value is arbitrary. There is some, but by 
no means universal, consensus that the effect of missing higher-order terms can 
hence be estimated from the dependence of as(Af|-) on the value of fi assumed 
in fits of the calculations to the data, and a renormalisation scale uncertainty is 
sometimes quoted. This procedure, well-motivated in that the //-dependence caused 
by the truncation of the perturbation series would be cancelled by addition of the 
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higher-order terms, is, however, arbitrary, and is not equivalent to knowledge of 
the size of the a priori unknown terms. In cases where scale uncertainties are 
considered this arbitrariness is manifested in the wide variation among the ranges 
and central values of // chosen by different experimental groups, see eg. [8]; in other 
cases this source of uncertainty is not included in the errors. Different as{M'^) results 
with similar experimental precision can hence be quoted with different total errors 
depending on the procedure adopted for assigning the theoretical uncertainties. The 
interpretation of the central values and errors on as(Af|-) measurements is hence not 
always straightforward. In this review it is attempted to summarise the sources of 
error and assumptions made in the various determinations. 

3 e+e~ Annihilation 

Since hadronic activity is restricted to the final state, the simplest environment 
for measurement of is provided by e"'"e~ annihilation. Experimental signatures 
of hadronic events are largely free of backgrounds, and the smaller number of 
Feynman diagrams contributing at a given order in perturbation theory renders QCD 
calculations more tractable than in lepton-hadron or hadron-hadron collisions. 

3.1 R and Lineshape 

For the inclusive ratio R = cr[e~^e~ — )• hadrons)/(T(e~'~e~ — )• the SM electroweak 

contributions are well understood theoretically and the perturbative QCD series has 
been calculated up to O(a^) [9] for massless quarks and up to O(a^) including quark 
mass effects [10]; the large size of the O(a^) term is potentially a cause for concern 
about the degree of convergence of the series. The SM prediction for R with the top 
quark mass fixed to 179 GeV, Mz fixed to 91.185 GeV, and the SM Higgs mass Mh 
fixed to 300 GeV was fitted to published measurements below the Z° resonance, in 
the cm. energy range 5 < Q < 65 GeV, and yielded [II]: 

a,(M|) = 0.mtlf,l{exp.)±0m2{MH). (4) 

The experimental error is large due to the small data samples and/or limited precision 
of the luminosity measurements. The second error has been assigned for variation of 
Mh between 60 and 1000 GeV. 



4 



The CLEO Collaboration is currently attempting to measure R at Q = 10.5 GeV. 
They expect to achieve a precision of order AR ~ ±2.9% [12], which will translate 
into a precision of about ±0.04 on as{M'^)^ providing a useful, but not very precise, 
low-energy measurement. 

Closely-related observables at the Z° resonance are the Z° total width F^, the 
pole cross section (7°, and the ratio of leptonic to hadronic Z° decay branching ratios 
i?;, which have been measured using the sample of approximately 16M hadronic 
Z° decays measured at LEP. In these cases the non-perturbative contributions 
are expected to be 0(1/ Mz) and are usually ignored. A concern is that recent 
measurements of observables that probe the electroweak couplings of the Z° to b 
and c quarks deviate slightly from SM expectations [13]. Since these couplings must 
be known in order to extract as(Af|-), this effect, whatever its origin, is a potential 
source of bias [3]. Eurther analysis is in progress from the SLC and LEP experiments 
and the situation is not yet resolved. 

Proceeding nonetheless, the procedure adopted [13] is to perform a global fit to 
a panoply of electroweak data that includes the W and top quark masses as well 
as the Z° observables relating to the lineshape, left-right production asymmetry, 
decay fermion forward-backward asymmetries, branching ratios to heavy quarks, and 
T polarisation, by allowing Mh and as{M'^) to vary. Data presented at the 1996 
summer conferences yield the positively-correlated results [13] Mh = 149^82° GeV 
and 



The as{M'^) value is lower than the corresponding results presented at the 1995 
conferences [14], a,(M|) = 0.123±0.005, and at the 1994 conferences, a,(M|) 0.125± 
0.005 [15], whose large central values were partly responsible for a supposed 
discrepancy between 'low-Q^' and 'high-Q^' ag^M^) measurements [16]. The change 
between 1995 and 1996 is due to a combination of shifts in the values of the Z° 
lineshape parameters, redetermined in light of the recalibration of the LEP beam 
energy due to the 'TGV effect' [13], and a change in the central value of Mh at 
which a,(M|) is quoted, from 300 GeV (1995) to the fitted value 149 GeV (1996). A 
detailed study of theoretical uncertainties implies [17] that they contribute at a level 
substantially below ±0.001. Since data-taking at the Z° resonance has now been 
completed at the LEP collider the precision of this result is not expected to improve 
further. 




(5) 
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3.2 r Decays 

An inclusive quantity similar to R is the ratio i?^ of liadronic to leptonic decay 
branching ratios, and Bi respectively, of the r lepton: 



where Be and B^ can either be measured directly, or deduced from a measurement of 
the T lifetime r^. In addition, a family of observables known as 'spectral moments' of 
the invariant mass-squared s of the hadronic system has been proposed [18]: 



where Af^ is the r mass. In this case the integrand can be measured independently 
of Be- It is easily seen that = i?™. 

i?^ and R'^' have been calculated perturbatively up to O(a^). However, because 
M-r ~ 1 GeV one expects as(Af^) ~ 0.3 and it is not a priori obvious that the 
perturbative calculation can be expected to be reliable, or that the non-perturbative 
contributions will be small. In recent years a large theoretical effort has been devoted 
to this subject; see eg. [18, 19, 20]. 

The ALEPH Collaboration derived i?^ from its own measurements of i?e, -B^, and 
T^, and also measured the (10), (11), (12), and (13) spectral moments. A combined 
fit yielded [21] a,{M^) = 0.330 ± 0.046 (a.(M|) = 0.118 ± 0.005) and = 
0.003 ±0.005. A later conference report [22] quoted a, (M|) = 0.124 ± 0.0022 ± 0.001, 
where the first error receives equal contributions from experiment and theory, and 
the second derives from uncertainties in evolving across the c and b thresholds; 
this result will be included in the average presented here. The OPAL Collaboration 
measured i?^ from Be, i?^, and r^, and derived [23] as(Af^) = 0.375to!oi8 (exp.) to!oi7 
(pert.) ±0.006 (non-pert.) (a.(M|) =0.1229tEJ:EJEJl? (exp.) f^omll (theor.)). The 
CLEO Collaboration measured the same four spectral moments as ALEPH and also 
derived i?^ using 1994 Particle Data Group values for i?e, B^ and r^. A combined fit 
yielded [24] a,(M^) = 0.306±0.017 (exp.) ±0.017 (theor.) (a.(M|) = 0.114±0.003). 
This central value is slightly lower than the ALEPH and OPAL values. If more recent 
world average values of Be and B^ are used CLEO obtains a central as(Af^) value 
of 0.339 [24]. Averaging the second CLEO result and the ALEPH and OPAL results 
by weighting with the experimental errors, assuming they are uncorrelated, yields 



Rr 



Bi 



1-Be-B, 

Be 



(6) 
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as{M'^) = 0.122 lb 0.001 (exp.) ±0.002 (theor.). This is nominally a very precise 
measurement. 

Recently, however, theoretical studies have been performed which suggest that 
additional uncertainties need to be added to the above result before it can be 
compared with other measurements of as{M^). A comparison of the exact 
O(a^) perturbative calculation of the spectral moments with a resummed procedure 
and with a partial resummation based on 'renormalon chains' concluded [20] that 
the theoretical uncertainty on as(Af^) is ±0.05. A study of ultraviolet renormalon 
ambiguities in the behaviour of the perturbative series yielded [25] an estimate of 
the theoretical uncertainty on as(Af^) of at least ±0.06. Until there is a resolution 
of these theoretical considerations it seems appropriate to conclude that the current 
value of as(Af|-) from r decays is: 

a,{Ml) = 0.122 ± 0.001(exp.) ± 0.006(theor.) (8) 
3.3 Hadronic Event Shape Observables 

In e"'"e~ annihilation one can define infra-red- and coUinear-safe measures of the 
topology of hadronic final states. Perhaps the most familiar such quantity is the rate 
of production of 3-jet events, defined using an iterative particle clustering algorithm, 
see eg. [26], but there are other measures relating to longitudinal momentum fiow in 
events, jet masses, and energy-energy correlations between particles; for a discussion 
see eg. [27]. The observables are constructed to be directly proportional to at 
leading order, and so are potentially sensitive measures of the strong coupling. To 
date the leading and next-to-leading order terms have been calculated for a large 
number of observables [28] . 

The technology of this approach has been developed over the past 15 years 
of analysis at the PETRA, PEP, TRISTAN, SLC and LEP coUiders, so that the 
method is considered to be well understood both experimentally and theoretically. 
It is necessary to correct the measured distributions for any bias effects originating 
from the detector acceptance, resolution, and inefficiency, as well as for the effects 
of initial-state radiation and hadronisation, to yield 'parton-level' distributions 
which can be compared directly with the QCD calculations. In the absence of 
non-perturbative calculations of the hadronisation process, Monte Carlo models of 
jet fragmentation [29, 30] are commonly used for this purpose. 

The Z° experiments have measured as(Af|-) in this fashion; see [8] and 
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references therein. For example, the recent SLD study employed 15 observables 
and found considerable scatter among the 15 as(Af|-) values, the r.m.s. deviation 
being 0.0076. This scatter is much larger than the hadronisation uncertainties 
and can be interpreted as arising from the effects of the uncalculated higher-order 
perturbative contributions. In fact, within the renormalisation scale uncertainties 
the as(Af|-) values from all the observables are consistent. The experimental 
collaborations have chosen different sets of obervables and different ranges over which 
to vary in combination with their different averaging methods this has led to some 
variation among the quoted central values of as(Af|-) and scale uncertainties [8]. 

For six observables improved calculations can be formulated that incorporate the 
resummation [31] of leading and next-to-leading logarithmic terms matched to the 
0{al) results. The matched calculations are expected a priori both to describe the 
data in a larger region of phase space than the fixed-order results, and to yield a 
reduced dependence of on the renormalization scale; they have been applied by 
the Z° experiments to determine as(Af|-) [27, 32]. 

Hinchliffe has reviewed the various measurements from experiments performed 
in the cm. energy range 10 < Q < 91 GeV, utlilising both O(a^) and resummed 
calculations, and quotes an average value of as(Af|-) = 0.122 it 0.007 [3], where the 
large error is dominated by the scale uncertainty, which far exceeds the experimental 
error of about ±0.002. Schmelling has also compiled the measurements, including the 
recent results from the LFP run at Q ~ 133 GeV [33], and quotes a global average 
[34] as(Af|-) = 0.121 it 0.005, in agreement with [3], but assuming a more aggressive 
scale uncertainty. The value 

a,(M|) = 0.121 ± 0.002(exp.) ± 0.005(theor.) (9) 

will be used here. 

The best way to reduce the theoretical uncertainty would be to calculate the 
observables to higher order in perturbation theory, a difficult task that has not yet 
been performed. In the absence of O(a^) QCD calculations it has been suggested 
[35] that the O(a^) calculation for each observable can be 'optimised' by choosing a 
specific value of the renormalisation scale. This approach has recently been applied 
[8] to 15 event shape observables using the 'PMS', 'FAC, and 'BLM' optimised 
scales. However, for any of these scale choices the scatter among the as(Af|-) values 
is comparable with that from the choice // = implying that higher-order effects 
contribute roughly equally in all of these procedures. Recently Fade Approximants 
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have been applied to estimate the O(a^) contributions to the QCD series for the same 
15 event shape observables [36]. Remarkably the scatter among the a^^M^) values is 
noticeably smaller than in the O(a^) cases, the r.m.s. deviation being ±0.0035. This 
convergence implies that the Pade method provides at least a partial approximation 
of higher-order perturbative QCD contributions, but explicit calculation of the 
O(a^) terms will be necessary to confirm this. 

Finally, the L3 Collaboration has utilised events with a hard radiated final-state 
photon, which reduces the effective cm. energy available to the hadronic system, 
to examine the Q^-evolution of four event shape observables in the range 30 < 
Q < GeV. By comparing with resummed + O(a^) calculations they derived [37] 
preliminary values of in each energy bin that, when evolved to Af^, are consistent 
with a,(M|) = 0.120. 

3.4 Scaling Violations in Fragmentation Functions 

Though distributions of final-state hadrons are not, in general, calculable in 
perturbative QCD, the Q^-evolution of the scaled energy [xp = 2E/Q) distributions 
of hadrons, or 'fragmentation functions', can be calculated and used to determine 
as- In addition to the usual renormalisation scale a factorisation scale fip must 
be defined that delineates the boundary between the calculable perturbative, and 
incalculable non-perturbative, domains. Additional complications arise from the 
changing composition of the underlying event fiavour with Q due to the different 
Q-dependence of the 7 and Z° exchange processes. Since B and C hadrons typically 
carry a large fraction of the beam momentum, and contribute a large multiplicity 
from their decays, it is necessary to consider the scaling violations separately in b, c, 
and light quark events, as well as in gluon jet fragmentation. 

In an early analysis [38] the DELPHI Collaboration parametrised the 
fragmentation functions using the O(a^) matrix elements and the string 
fragmentation model implemented in JETSET [29]. They fitted data in the range 
14 < Q < 91 GeV to determine ag^M^) = 0.118ib0.005, where the error is dominated 
by varying fi in the range 0.1 < fi/Q < 1. The ALEPH Collaboration used its 
Z° data to constrain fiavour-dependent effects by tagging event samples enriched 
in light, c, and b quarks, as well as a sample of gluon jets [39]. The fragmentation 
functions for the different fiavours and the gluon were parametrised at a reference 
energy, evolved with Q according to the perturbative DGLAP formalism calculated 
at next-to-leading order [40], in conjunction with a parametrisation proportional 
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to IjQ to represent non-perturbative effects [41], and fitted to data in the range 
22 < g < 91 GeV. Tliey derived a,(M|) = 0.126 ± 0.007 (exp.) ±0.006 (tlieor.), 
where the theoreticai uncertainty is dominated by variation of the factorisation scaie 
IJ.p in the range — 1 < Infi^/Q^ < 1; variation of the renormalisation scale in the same 
range contributed only ±0.002. DELPHI has recently reported a similar analysis 
[42] yielding a,(M|) = 0.12lto:oo7 (exp.) ±0.010 (theor.). Curiously, although a 
similar range as ALEPH, 0.3 < fi/Q < 3, was used to examine variation of the 
renormalisation and factorisation scales, here the renormalisation scale dominates 
the theoretical uncertainty, with a contribution of ±0.009, in contrast to ±0.002 
from factorisation. Combining the ALEPH and later DELPHI results, assuming 
uncorrelated experimental errors, yields: 

a,{Ml) = 0.124 ± 0.005(exp.) ± 0.010(theor.) (10) 

where, until the apparent discrepancy between the ALEPH and DELPHI analyses of 
theoretical uncertainties is resolved, the larger DELPHI value is taken. 

4 Lepton-hadron Scattering 

The study of scaling violations in structure functions has been of historical significance 
in establishing QCD as the theory of strong interactions. Measurements have been 
performed using electron, muon and neutrino beams directed onto a variety of 
both polarised and unpolarised nuclear targets. An enormous range of the (x, Q^) 
kinematic plane has thus been probed, where Q is the four-momentum transfer and 
X the fractional proton energy carried by the struck parton. This continues to be an 
active field, with the HERA electron-proton collider at DESY producing interesting 
results at the lowest-x and highest-Q^ values yet reached, and providing new tests of 
QCD from the jet structure of the hadronic current. 

An important theoretical issue is the presence of power-law corrections to the 
perturbative Q^-evolution, or 'higher-twist' contributions, of the form which 
enter with a priori unknown x-dependent coefficients that cannot in general be 
calculated. In cases where higher-twist contributions are considered, the coefficient 
is usually extracted by fitting to the data and is correlated with the value thus 
determined. A corresponding uncertainty on can then be assigned from its 
dependence on the size of the higher-twist contribution varied within some reasonable 
range allowed by the data. 
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4.1 Unpolarised Structure Functions 

The nucleoli structure function can be decomposed into pieces labelled F2{x^Q^) 
and xF3[x^Q^)^ each of which may be written as a sum of singlet and non-singlet 
contributions with respect to the quark flavour content, the former having no net 
flavour. Except at low-x, only the evolution, not the actual shape, of the structure 
functions is predicted by QCD, and it is different for the singlet and non-singlet 
components. In the non-singlet case the evolution enters only via as{Q^)^ whereas 
in the singlet case it enters via both as{Q^) and, for x < 0.5, the gluon distribution 
function g{x, Q^). The dependence of the F2 structure function has been measured 
in electron, muon, neutrino and antineutrino scattering experiments. The singlet 
contribution has been determined using deuterium (e, //) and iron (p) targets, and 
the non-singlet contribution using hydrogen (e, //) targets. The dependence of 
the non-singlet XF3 structure function has been measured using neutrino scattering 
on iron. 

A detailed analysis of the Q^-evolution of measurements of F2 was presented 
several years ago [43] , in which a global fit was performed to hydrogen and deuterium 
data from electron scattering at SLAC and muon scattering at CERN. Eor the 
'non-singlet' region, x > 0.25, significant higher-twist contributions were required to 
fit the data. The effect of these terms was reduced by restricting the fit to high-Q^ 
bins to yield: 

a,{Ml) = 0.113 ± 0.003(exp.) ± 0.004(theor.) (11) 

where the 'experimental error' includes normalisation uncertainties on the data, and 
uncertainties on the higher-twist and gluon distribution contributions; the theoretical 
uncertainty is dominated by varying the renormalisation and factorisation scales in 
the range 0.1 < fi'^/Q^ < 4. 

The NMC Collaboration also determined using their F2 data, but at low-x 
[44]. In this case the evolution is dominated by the gluon distribution, which 
must be constrained simultaneously with a^. Higher-twist terms were included as 
fixed contributions and were not determined from the data. NMC found: as{7 GeV^) 
= 0.264 ± 0.018 (stat.) ±0.070 (syst.) ±0.013 (higher-twist); no uncertainty due to 
scale variation was assigned. This translates to: 

a,{Ml) = 0. 117 ± 0.003(stat.)+°:°J°(syst.)± 0.002(higher twist). (12) 

The CCER Collaboration similarly fitted the Q^-dependence of XF3. They 
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required >15 GeV^ in order to minimise higher-twist effects, and found [45]: 
a,{Ml) = 0.111 ±0.002 (stat.) ±0.003 (syst.) ±0.004 (theor.), where the systematic 
error is iargeiy due to the uncertainty in the energy caiibration of the detector. This 
result has very recently been updated following a recalibration of the energy scale. A 
simultaneous fit to F2 and XF3 yields [46]: 

a,{Ml) = 0.119 ± 0.0015(stat.) ± 0.0035(syst.) ± 0.004(theor.); (13) 

it can be seen that the central value has increased by 0.008 relative to the published 
result. 

CCFR has also studied the Q^-dependence of the integral over x of XF3 that defines 
the Gross-Llewellyn-Smith sum rule [46]. This quantity is suitable for determination 
of as as it is independent of the gluon distribution, and it has been calculated up to 
O(a^) in perturbative QCD. The preliminary CCFR measurement is: 

a,(M|) = 0.108i°:°°^(stat.) ± 0.004(syst.)l°:°°;^(higher twist). (14) 

An update of this result that refiects the energy recalibration of the detector is eagerly 
awaited. 

Finally, has recently been determined in a theoretical study of the F2 structure 
function of the proton, using HFRA data at low x and large [47]. In this kinematic 
region the structure function exhibits 'double scaling', suggesting that a resummation 
of terms containing powers of logQ^ and logl/x can be compared reliably with the 
data. Such a resummed calculation has been performed up to next-to-leading order 
in as including all leading and subleading logarithmic terms. Higher-twist corrections 
are also claimed to be small. A comparison with the 1993 HFRA data yields [47]: 

a,(M|) = 0.120 ± 0.005(exp.)± 0.009(theor.). (15) 

In this case the largest contribution to the 'experimental' error actually derives 
from uncertainties associated with the parton distribution functions (PDF), and 
the theoretical uncertainty is dominated by changing the renormalisation and 
factorisation scales in the range 0.5 < fi/Q < 2. Since this technique appears to 
offer the promise of a reasonably precise determination of that is independent of 
other methods used in deep-inelastic scattering, one hopes that it will be adopted 
by the HFRA experimental collaborations and applied to their own data so that 
experimental errors can be properly assessed. 
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4.2 Polarised Structure Functions 



In deep-inelastic scattering ol polarised leptons on polarised nuclear targets the 
longitudinal spin structure functions of the proton and neutron can be measured. 
The Bjorken sum rule is defined by the integral over x of the difference between these 
structure functions; from the point-of-view of the determination of this quantity is 
similar to the GLS sum rule and it has also been calculated perturbatively up to O(a^). 
Data from a variety of experiments have been examined by the E143 Collaboration 
[48] and the Q^-dependence of the sum rule has been used to determine: as(Af|-) = 
0.119lo!oi9 (exp.) assuming no higher-twist contribution, and as(Af|-) = 0.113lo!o35 
(exp.) when an ad hoc parametrisation of such a term was included in the fit; no 
additional theoretical uncertainties were considered. 

A subset of these data has been included in detailed studies of higher-order 
perturbative QCD contributions to the Bjorken sum rule. Including estimated 
O(a^) contributions, and neglecting higher-twist effects, the value as(Af|-) = 
0.122lQ;QQg (exp.) was obtained [49]. Using a QCD sum-rule estimate of the 
higher-twist coefficient, with a ±50% uncertainty, and assigning an uncertainty 
of ±10% to the O(a^) coefficient, adds an additional theoretical uncertainty and 
leads to a,(M|) = O.llStoml 

[49]. A further recent analysis, including a study 
of the renormalon structure of the Borel transform of the perturbative prediction 
for the Bjorken sum rule, as well as application of Pade approximants to estimate 
O(a^) contributions, yielded [50]: 

a,{Ml) = 0.116l°:°°^(exp.) ± 0.003(theor.) (16) 

where the theoretical uncertainty includes contributions from the uncertainty on the 
Pade estimate, as well as from the higher-twist coefficient uncertainty. This nominally 
precise result is encouraging, and warrants application of these new theoretical 
techniques by the experimental collaborations to their own data. 

4.3 Jet Final States 

Hadronic final states in lepton-hadron scattering can be analysed in a similar 
fashion to those in e"'"e~ annihilation (Section 3.3) via infra-red- and coUinear-safe 
measures. The ratio of cross sections i?2-i-i = C2-1-1 /ctot, where the '2' refers to events 
containing two resolved jets in the final state in addition to '±1' proton remnant 
jet, is particularly suitable since it is directly proportional to and has recently 
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been calculated at next-to-leading order [51] for the JADE jet algorithm [52]. This 
observable is subject to similar hadronisation uncertainties as those encountered in 
e"'"e~ annihilation, and to additional sources of uncertainty relating to the parton 
distribution functions of the proton as well as contributions to hadronic activity from 
initial-state parton interactions. 

The HI Collaboration divided their data sample into 5 bins of in the range 
10 < < 4000 GeV^ and evaluated i?2+i. Using PDFs fitted to HERA data, they 
compared the QCD prediction with data in the two highest bins and measured [53]: 
a,(M|) = 0.123 ± 0.012 (stat.) ±0.007 (syst.) ±0.011 (theor.). The experimental 
systematic error is dominated by the uncertainty on the hadronic energy scale of the 
calorimeter, and the theoretical uncertainty includes contributions from the QCD 
model-dependence of the acceptance corrections, choice of structure functions, value 
of j/c, and variation of the renormalisation and factorisation scales in the ranges 
1/4 < fJ'^ I < 4. A similar analysis by the ZEUS Collaboration, based on 3 
bins in the range 120 < < 3600 GeV^, found [54]: a,(M|) = 0.117 ± 0.005 
(stat.) toZt (syst.) toZl (had.) ±0.001 (PDEs) toZl (scale), where variation of the 
renormalisation and factorisation scales in the range 0.4 < fi^ /Q^ < 2 was considered. 
Averaging over these values assuming uncorrelated experimental errors yields: 

a,{Ml) = 0.118 ± 0.006(exp.)±0.011(theor.). (17) 

Interestingly ZEUS and HI differ in their estimates of theoretical uncertainties from 
nominally the same sources. Eor example, from PDEs ZEUS quotes ±0.001, whereas 
HI quotes ±0.005, and from variation of the renormalisation and factorisation scales 
ZEUS quotes to!oo6 HI only ±0.003, despite the fact that HI considered a larger 
range of scale change. Until this situation is resolved the larger HI estimate of the 
theoretical uncertainty has been applied to the average value. 

Since these are the first such analyses from HERA, and since similar analysis 
technology has been in use in e"'"e~ annihilation for over 15 years, there is every 
reason to expect that the sophistication of the analyses will improve and that the 
error estimates will stabilise. It will be interesting to see if different HERA jet 
observables are as sensitive to uncalculated higher-order perturbative contributions 
as those in e"'"e~ annihilation, and to what extent they yield a similar degree of scatter 
in a,(M|) (Section 3.3). 
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5 Hadron-hadron Collisions 



Hadron-hadron collisions provide a challenging environment for precise measurement 
of as- On the experimental side final-state jet measurements are complicated 
by remnants of the initial-state hadrons, and on the theoretical side calculations 
at next-to-leading order in perturbation theory include contributions from a large 
number of Feynman diagrams and have only recently been achieved for jet processes. 
Finally, the coefficients of powers of in the jet matrix elements depend upon 
the parton distribution functions of the incoming hadrons; these are constrained at 
low scales by deep-inelastic scattering experiments, but uncertainties accrue in 
their extrapolation to the higher jet transverse energy, Et^ scales that are currently 
accessible at the TeVatron. Any attempt to extract from the jet data is hence 
affected by these uncertainties, especially in the gluon distribution, and care must be 
taken to extract simultaneously a common value of used in both the hard-scattering 
matrix elements and the evolution of the PDFs. 

5.1 {W + l-iet)/{W + 0-jet) Ratio in pp Collisions 

Several experimental collaborations have attempted to determine from the ratio 
R of the cross sections for production of final states containing a W boson -|- 1-jet 
and a W boson -|- 0-jets; R is proportional to at leading order, and many sources 
of experimental uncertainty are expected to cancel. The first two studies of this 
technique were performed at the CFRN pp collider before a full next-to-leading 
order calculation of R had been achieved, and hence relied upon calculations of 
so-called 'A^-factors' to estimate higher-order QCD corrections to the leading-order 
result. Data from the UAl Collaboration were used to measure [55]: as{M^) 
= 0.127 lb 0.026 (stat.) ±0.034 (syst.), where the systematic error is dominated 
by a contribution of ±0.025 from the A^-factor uncertainty, with additional large 
contributions of ±0.013 from PDFs and ±0.010 from fragmentation modelling. The 
UA2 Collaboration performed a similar analysis, but used next-to-leading order 
calculations of the total W production cross section as well as the W px-distribution, 
to make a more accurate determination of the A^-factors. They found [56] as{M^) 
= 0.123 ± 0.018 (stat.) ±0.017 (syst.); the systematic error comprises components 
of ±0.012 from experimental sources and ±0.011 from theoretical uncertainties, the 
dominant contribution to the latter, ±0.010, arising from fragmentation modelling. A 
strong sensitivity to the renormalisation scale was found; variation of fi from Mw to 
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Mw I'i caused a change Aa^ = —0.010, but no corresponding theoretical uncertainty 
was assigned. 

Averaging these results by weighting by the experimental errors assuming they 
are uncorrelated, and evolving to the scale Af^, yields: 

a,{Ml) = 0.121 ±0.018(exp.) ±0.011(theor.) ±0.010(scale) (18) 

where the UA2 theoretical uncertainty has been taken and the last error has 
been assigned, based on the UA2 scale-dependence, as an estimate of higher-order 
uncertainties. 

Recently a complete next-to-leading order calculation of R has been performed 
[57], and has been utilised by the DO Collaboration in an attempt to measure [58]. 
Their measurement of R with an accuracy of about ±10% allows in principle a 
determination of as{M^) with an experimental error of ±0.016. However, they were 
unable to obtain a value of from a fit of the hard matrix element that is consistent 
with the value used as input to the parton distribution functions, and concluded that 
the sensitivity of this technique to is much smaller than expected. 

5.2 Inclusive Jet Cross Sections 

In principle can be determined from inclusive jet cross section measurements, for 
example the inclusive 1-jet Et distribution. A recent CDF measurement of the latter 
observable [59] has aroused much interest due to speculation on the origin of the 
population of events at high Et- This population is now believed to be explicable by a 
larger gluonic content of the proton at high fractional momentum than was previously 
expected [60]; this provides graphic illustration of the necessity to constrain parton 
distribution functions before other parameters, such as a^, can be determined from 
the data. 

A demonstration of the possibility to determine from the single-jet inclusive 
Et distribution has been provided recently [61]. Using a given value as input to a 
particular parametrisation of PDFs, and fitting the next-to-leading order prediction 
to the CDF Run la data for Et > 30 GeV, it appears possible to measure as{M'^) to 
a precision of ±0.001 (stat.) ±0.008 (syst.) ±0.005 (theor.), where the theoretical 
uncertainty was defined by variation of the renormalisation scale in the range 0.5 < 
fJ'l Et < 2, and potential non-perturbative effects were not considered. This is an 
encouraging first step. A programme of further studies is planned, whereby the 
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dependence on the choice of PDFs will be systematicUy studied, and ultimately the 
PDFs themselves will be extracted simultaneously with a^^M^) by fitting to the 
triply-difi'erential di-jet inclusive distributions [61]. 

Recently the UAl Collaboration has released a measurement of based on the 
cross section for bb production. They performed a fit of the next-to-leading order 
calculation of 'quasi-2-body' bb final-states to their data [62] and found: as(20 GeV) 
= 0.145tEJ:EJl^ (exp.) ±0.003 (m^) (scale) t^fw (PDF), which can be evolved to: 

a,{Ml) = 0.113i°:r(exp.)i°Z(theor.), (19) 

where the (dominant) renormalisation scale uncertainty corresponds to a variation 

1/4 < nlQ < 1 with Q = \J {krubY + pi and k is an additional fit parameter. 
This result represents one of the most precise determination of as{M'^) in hadronic 
collisions to date. 

5.3 Direct Photon Production in pp and pp Collisions 

Another reasonably precise determination was provided by the UA6 Collaboration, 
who measured inclusive direct photon production cross sections in pp and pp collisions 
using a molecular hydrogen gas-jet target at the CFRN pp collider [63]. They 
were able to isolate the contribution from the process qq — )• 7 g, which has been 
calculated at next-to-leading order in [64]. Using parton distribution function 
parametrisations fitted to BCDMS data they determined: 

a,{Ml) = 0.112 ± 0.006(stat.)± 0.005(syst.)+°:°°5'(theor.), (20) 

where the theoretical uncertainty was defined by variation of the renormalisation scale 
in the range p^/S < /U^ < 3pj/4, the central value being based on the PMS-optimised 
scale (Section 3.3). 

6 Heavy Quarkonium Systems 

Heavy quarkonium systems can be used to determine either from the measured 
hadronic decay rate or via the strength of the binding provided by the strong potential; 
the latter is achieved in practice by comparing the measured energy-level splittings 
with a lattice QCD calculation. 
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6.1 Heavy Quarkonium Decays 

The partial widths Tggg and T^gg lor the decay ol the states ol the J/^ and 
T into 3 gluons, or a direct photon + 2 gluons, respectively, have been calulated 
perturbatively at next-to-leading order. By dividing by the partial width Tn lor the 
leptonic decay the dependence on the wavelunction cancels and can in principle 
be determined. A non-relativistic analysis [65] revealed large differences between 
the values determined in this fashion from the J/^ and T systems, as well 

as a large renormalisation scale dependence. An improved analysis incorporated ad 
hoc relativistic corrections as well as next-to-next-to-leading order (NNLO) terms that 
were determined by fitting simultaneously to the data on Tggg/Tn for the J/^ and T 
IS - 3S resonances. The relativistic corrections turn out to be large for the J/^. A 
preliminary result: a, (Ml) = 0.113 ±0.001 ( exp.) ^Q^QQ5 (theor.) was obtained [65], 
where the theoretical uncertainty was defined from the variation of the fitted NNLO 
term for a scale variation in the range 0.5 < fi/Q < 2; a final version of this analysis 
has yet to be published. 

CLEO similarly used their measured ratio ^-ygg/^ggg for the T(1S) and have 
presented a prehminary result [12] as(M|) = 0.110 ±0.001 (exp.) ±0.004 (syst.) 
±0.005 (theor.), where the last uncertainty was defined by 0.1 < fi^/Q^ < 1. 
Combining with the analysis of [65], assuming uncorrelated experimental errors, yields 
a preliminary average: 

a,{Ml) = 0.113 ±0.001(exp.)+°:°°^(theor.). (21) 

6.2 Lattice Gauge Theory 

Lattice gauge theory currently provides a successful tool for performing non- 
perturbative QCD calculations, although it is presently limited in applicability to 
static properties of hadrons. Experimental data on hadron properties, such as meson 
masses or energy-level splittings in heavy mesons, can be input to lattice simulations 
and as extracted by calculation; for a recent review see [66]. Currently the precision 
of such determinations of as(Af|-) is not limited by experiment, but by the various 
theoretical uncertainties relating to lattice discretisation, treatment of 'sea quarks', 
and matching between the different renormalisation schemes used in lattice and 
perturbative calculations. 

Considerable interest was generated several years ago with the determination by 
the ENAL/SCRI group of a nominally very precise as(Af|-) measurement using the 
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IP - IS splitting in the charmonium system: as{M'^) = 0.105 it 0.004 [67]. The 
dominant error was derived from extrapolating from the 'quenched approximation', 
corresponding to zero quark flavours in vacuum polarisation contributions to the 
strong coupling, to four-flavour contributions above the charmonium energy scale. 
A somewhat larger, and even more precise, result, based on input data from the 
T spectrum, was presented more recently by the NRQCD group [68]: as(Af|-) = 
0.115 lb 0.002. In this case both zero and two flavours of dynamical quarks were 
considered, and the result extrapolated to three flavours, before evolving through 
the cc and bb mass thresholds up to Mz- The dominant error was derived from 
the uncertainty in the two-loop matching between the ay and MS schemes. An 
independent determination based on the p mass and the charmonium IS-IP splitting, 
and using two flavours of dynamical quarks, yielded [69]: as{M'^) = 0.108 it 0.008. A 
review of these measurements quoted [66]: as{M'^) = 0.112 it 0.007, the conservative 
error estimate reflecting a dispersion among the values from the different groups that 
is larger than the quoted errors in some cases. 

Further progress has been made recently; the two-loop matching coefficient has 
been calculated for the quenched approximation, and errors relating to the dynamical 
quark mass value have been considered. New preliminary as(Af|-) results have been 
presented [70] by the FNAL/SCRI group: 0.116 ± 0.003, and by the NRQCD group: 
0.118 lb 0.003, where in both cases the two-loop matching uncertainty is dominant. 
This recent convergence of the lattice results is encouraging, although the shifts in the 
central values of as(Af|-), as well as changes in the error estimates, as the techniques 
have improved is perhaps an indication that the systematics are not yet fully under 
control. An average over the two recent preliminary results yields: 



where the error is largely from common theoretical uncertainties. 

7 Summary and Outlook 

The average as(Af|-) value from each method, derived as described in the text, 
is shown in Table 1, together with the total experimental error and theoretical 
uncertainty. For the benefit of the aesthetically-minded the same data are summarised 
in Figure 1. Perhaps the most difficult aspect of this review is to combine these 
17 results in a meaningful way. It is worth restating that, in all cases, except 
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perhaps the Z° hneshape measurement, the theoretical uncertainty is a non-negligible 
contribution to the total error, and in many cases it is larger than the experimental 
error. By definition such theoretical uncertainties can only be estimated in an ad 
hoc fashion, and the degree of conservatism of the assigned uncertainty can often 
be contested. I have attempted to summarise the contributions to the theoretical 
uncertainty as quoted by the authors of the respective measurements and to assess 
the 'degree of reasonableness' of their estimate. In those cases where subsequent 
developments suggest a different estimate than that originally assigned, I have 
changed the theoretical uncertainty accordingly; in all cases the original results are 
included in the text and can be recovered and treated differently by those who hold 
a different view. 

It is certainly not correct to assume that all 17 measurements of 
as(Af|-) summarised in Table 1 are completely independent of one another. 
Correlations between measurements and between the systematic errors, especially 
on the theoretical side, certainly exist. However, no attempt has previously been 
made in a review of this kind to evaluate the correlations between the different 
measurements and to treat the errors accordingly in forming an average; this tradition 
will not be broken here. Taking an average over all 17 measurements assuming they 
are independent, by weighting each by its total error, yields as{M'^) = 0.118 with 
a of 6.4; the low value refiects the fact that most of the measurements are 
systematics-limited. Taking an unweighted average, which in some sense corresponds 
to the assumption that all 17 measurements are completely correlated, yields the same 
result. There can hence, hopefully, be little dispute that the average value as(Af|-) = 
0.118 accurately characterises the 17 measurements. 

The first procedure yields an error on the average of ±0.001; because of 
correlations this is almost certainly an underestimate of the 'true error'. In the spirit 
of the second procedure one could take the error on the most precise measurement, 
namely ±0.003; this has been adopted by a previous reviewer [3]. Instead I choose 
to quote the r.m.s. deviation of the 17 measurements, w.r.t. their average value, to 
characterise the dispersion: 



which is presumably a conservative estimate of the error. In an independent review 
[34] an attempt was made to consider positive correlations between measurements by 
rescaling errors according to the respective value. This yielded the same average 
as(Af|-) value as presented here, and an error of ±0.003. 




0.118 ± 0.005, 



(23) 
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It is self-evident from Fig. 1 that, within the errors, there is no evidence of 
any discrepancy between measurements made at 'low'- and 'high'-Q^. It is notable, 
however, that the precision of the individual measurements varies between about 3 
and 20%, so that any anomalous effects up to the 5% level would not necessarily be 
discerned. 

If further progress is to be made in testing QCD, future measurements of 
as(Af|-) should aim for substantially improved precision. The prospects for achieving 
1%-level measurements are discussed elsewhere [71]. Lattice QCD determinations 
may reach this precision within the next few years. Deep-inelastic scattering and 
e"'"e~ annihilation will probably require higher-energy facilities, as well as significant 
theoretical effort in the areas of higher-twist and O(a^) perturbative contributions, 
respectively. A precise as(Af|-) measurement has yet to emerge from the TeVatron, 
but feasibility studies are in progress and appear promising. 
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Method 


Q (GeV) 


a,(M|) 


exp. 


theor. 


e+e" : r decays 


1.8 


0.122 


0.001 


0.006 


DIS: Bjorken SR 


1 - 4 


0.116 


-|-0 . 003 
-0.005 


0.003 


DIS: GLS SR 


1 - 5 


0.108 


_|_Q QQ5 

-0.006 


_|_Q QQ4 

-0.006 


DIS: F2 (NMC) 


1 - 7 


0.117 


+ 0.010 
-0.015 


0.002 


DIS: F2 (HERA) 


1 - 10 


0.120 


0.005 


0.009 


DIS: F2 (SLAG, BCDMS) 


2 - 16 


0.113 


0.003 


0.004 


pp, pp: direct 7 


4 - 6 


0.113 


0.008 


+0.009 
-0.001 


LGT: J/*, T 


3 - 11 


0.117 


— 


0.003 


J/^, T decays 


3 - 11 


0.113 


0.001 


+ 0.007 
-0.005 


DIS: F2, xFs (CCFR) 


3 - 23 


0.119 


0.004 


0.004 


DIS: jets 


10 - 63 


0.118 


0.006 


0.011 


pp: bb prod. 


20 


0.113 


+ 0.007 
-0.006 


+ 0.008 
-0.009 


e+e": R 


5-65 


0.128 


+ 0.012 
-0.013 


0.002 


pp: W+l-jet 


80 


0.121 


0.018 


0.015 


e+e": event shapes 


10 - 133 


0.121 


0.002 


0.005 


e+e": fragmentn. fns. 


14 - 91 


0.124 


0.005 


0.010 


e+e~: Z° lineshape 


91 


0.1202 


0.0033 


< 0.001 



Table 1: Average as(Af|-) values and errors from the 17 methods 
described in the text. The approximate range in Q is shown for each 
method. 
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Figure 1: Average as(Af|-) values and errors from the 17 methods 
described in the text. The results are ordered vertically in Q. 
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